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This  report  will  be  broken  into  two  parts.  The  first  part 
will  deal  with  the  theory  underlying  the  principles  of  wavelength 
independent  optical  lithography  and  microscopy  and  the  second 
part  will  deal  with  the  experimental  realization  of  these 
theoretical  principles. 
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abstract 


The  transmission  of  light  through  an  infinite  slit  in  a  thick,  perfectly 
conducting  screen  is  investigated.  The  spatial  distribution  of  the  near-field 
energy  flux  Is  determined  through  the  formulation  of  four  coupled  integral 
equations,  which  are  solved  numerically.  Transmission  coefficients  calculated 
by  this  method  are  in  agreement  with  those  determined  by  an  alternative 
formulation.  The  results  theoretically  demonstrate  the  feasibility  of 
near-field  super-resolution  microscopy,  in  which  the  collimated  radiation 
passed  by  an  aperture  is  used  to  circumvent  the  diffraction  limit  of 
conventional  optics,  and  further  suggest  the  feasibility  of  near-field 
super-resolution  acoustic  imaging. 

INTRODUCTION 

The  transmission  of  radiation  through  small  apertures  has  been 
investigated  for  a  number  of  years.  Bethe1  treated  tee  case  of  diffraction  cy 
a  round  hole  in  an  infinitely  thin,  perfectly  conducting  screen  in  the  limit 
where  the  hole  was  small  compared  to  the  wavelength,  and  considered  the 
applicability  of  his  solution  to  frequency  behavior  and  power  flow  in 
microwave  cavities.  More  recently,  research  in  aperture  diffraction  has  beer, 
aimed  at  obtaining  a  better  understanding  of  tr.e  effects  of  nucie.  r-induced 
electromagnetic  pulses.1  However,  these  studies  nave  not  dealt  with  the 
near-field  spatial  distribution  of  radiation  transmitted  past  the  aperture. 

The  development  by  Ash  and  'iicrolls1  of  a  scanning  microwave  microscope  with 
i/60  resolution  provides  tr.e  incentive  to  apply  aperture  diffraction  theory  to 
a  new  problem:  tr.e  ietermi.-at.cn  of  tr.e  practical  limits  of  near-field 
super-resolution  microscopy.  trough  in  tms  paper  we  will  consider  the 
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specific  case  of  optical  microscopy,  the  results  are  applicaole  to  other 
regimes  of  the  electromagnetic  spectrum. 

In  the  near-fleld  microscope,  an  illuminated  subwavelength  aperture  is 
used  as  a  light  source.  This  aperture  is  then  placed  within  the  near  field 
relative  to  an  object  and  an  image  is  generated  by  scanning  the  aperture 
relative  to  the  object  (see  Fig.  1).  From  a  theoretical  standpoint,  the  most 
interesting  questions  concern  the  degree  of  collimation  of  the  radiation 
immediately  past  the  screen  and  the  rate  of  spreading  of  the  radiation  farther 
from  the  screen,  which  would  determine  the  resolution  capabilities  and 
positioning  requirements  of  a  super-resolution  instrument.  To  address  these 
questions,  a  two-dimensional  aperture  model  will  be  considered,  l.e.,  an 
infinite  slit  in  a  thick  screen  of  perfect  conductivity. 

The  infinite  slit  geometry  has  already  been  considered  by  a  number  of 
researchers.  Harrington  and  Auckland*  have  analyzed  the  problem  using  the 
method  of  moments  and  have  developed  a  simple  equivalent  circuit  in  the  narrow 
slit  limit  (width  <  x).  However,  their  results  are  concerned  primarily  with 
transmission  coefficients  rather  than  the  spatial  distribution  of  the 
radiation.  Neerhoff  and  Mur*  have  used  a  Green's  function  formulation  for  the 
problem.  Indeed,  their  derivation  will  oe  used  for  the  analysis  here,  and  is 
summarized  below.  The  difference  between  this  work  and  that  of  Neerhoff  and 
Mur  lies  in  the  manner  in  which  the  model  is  applied.  First,  the  smallest 
slits  considered  by  Neerhoff  and  Mur  were  roughly  A / 2  in  width,  whereas  much 
smaller  slits  (width  »X/tg;  are  of  interest  in  super-resolution  applications. 
Second,  although  the  Luthers  compare  the  electromagnetic  field  amplitudes  in 
the  plane  of  the  slit  wit.n  similar  amplitudes  farther  from  the  screen,  the 
down-screen  distances  c.nosen  for  comparison  are  generally  beyond  the  range  of 
collimation  as  determined  by  the  more  thorough  calculations  included  in  this 
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work.  Finally,  when  considering  illumination  for  microscopic  or 
photolithographic  applications,  the  energy  flux  transmitted  by  the  slit  is  of 
greater  physical  significance  than  are  the  field  amplitudes  calculated  by 
Meerhoff  and  Mur. 

In  short,  the  theoretical  development  of  Meerhoff  and  Mur  will  oe  used  to 
find  the  fields  at  the  boundaries  of  a  number  of  narrow  slits  and  these  fields 
will  then  be  used  to  calculate  the  spatial  distribution  of  the  transmitted 
energy  flux  within  the  near  field.  To  check  the  consistency  of  the  results, 
transmission  coefficients  will  be  calculated  and -compared  to  the  findings  of 
Harrington  and  Auckland.  The  implications  of  the  results  for  super-resolution 
microscopy  will  then  be  discussed. 

THEORETICAL  DEVELOPMENT 

Initial  Formulation 

The  following  formulation  summarizes  the  work  of  Meerhoff  and  Mur.  The 
notation  used  has  also  been  largely  adopted  from  that  work. 

The  configuration  appropriate  for  the  formulation  of  the  problem  is  shown 
in  Fig.  2.  The  slit,  with  a  width  2a,  is  chosen  parallel  to  the  y  axis  and  is 
contained  in  a  screen  of  thickness  b.  The  space  surrounding  the  slit  is 
divided  into  three  regions.  Region  1  (  j  x  |  < • .  jy|<»,  b<zo)  contains  the  space 
above  the  sLit  and  includes  one  incident  electromagnetic  wave.  Region  II 
;]xj<a,  |yj<«,  3<z<b)  contains  the  space  inside  the  slit.  Finally,  region  III 
( |x| <-,  |y|<«,  -«<z<0)  contains  tne  3pace  below  the  slit  and  includes  the 
transmitted  wave.  Different  oie^ectrics  can  exist  m  tne  different  regions. 

The  incident  pla^e  -ave  m  region  I  propagates  in  the  xz  plane  at  an 
angle  3  with  -espect  to  tne  o  oxis.  The  magnetic  fieli  is  assumed  to  be  time 
r.armonic  and  botn  polarized  a.'d  constant  in  tne  y  direction: 

M(x,y,z,t)  -  ’J(x,z)exp  e..  {',) 
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These  restrictions  are  placed  for  two  reasons.  First,  elementary  waveguide 
calculations  Indicate  tnat  at  least  one  propagating  mode  exists  for  this 
polarization,  no  matter  how  narrow  the  slit.  Hence,  the  arrangement  is  of 
physical  interest  since  large  transmission  coefficients  can  be  anticipated. 
Second,  Maxwell's  equations  can  be  used  to  find  the  electric  field  from 
U(x,z) : 


E  (x,z) 

-  3  U(  x ,  z  ) 

(2a) 

X 

wC  z 

E  ( x ,  z ) 

-  0 

(2b) 

E  (x,z) 
z 

-  —  3  U  ( x ,  z ) 

41 C  X 

(2c) 

Thus,  the  restrictions  in  Eq.  (1)  reduce  the  diffraction  proDlem  to  one 
involving  a  single  scalar  field  in  only  two  dimensions. 

The  field  is  represented  by  Uj(x,z)  ( J -1 ,2,3)  in  each  of  the  three 
regions.  The  time  harmonic  behavior  indicates  that  toe  field  satisfies  the 
Helmholtz  equation: 

(7*-k* )U  -0  (J-1.2,3)  (3) 

j  j 

In  region  I  the  field  is  further  decomposed  into  three  components: 

IVx.z)  «  U^Cx.z)  *  Jr(x,z)  *  ’J  3  ( x ,  z )  ;3 

each  of  which  satisfy: 


o  7  r  d 

.7 


(A  is  represents  the  incident  field,  which,  for  simplicity,  is  assumed  to  ce 
a  plane  wave  of  unit  amplitude: 

.,i.  ,  Ik  ,  ( xs  -ccos9  )  , - 

J  \x,z)-e  J 

'Jr  denotes  the  field  that  ce  -“fleeted  if  tr.ere  -ere  no  slit  m  t.ne 


screen,  and  thus  satisfies-. 


0 


■Jr(x,z)  -  Ul(x,2b-z) 


Ud  describes  ttts  diffracted  field  in  region  I  due  to  the  presence  of  the 


slit. 

For  a  perfectly  conducting  screen,  the  tangential  S  vanishes  at  the 
surface,  so  that: 

3nU-0  at  the  conductor  surface  (9) 

where  n  is  the  normal  vector  at  the  surface.  Further  boundary  conditions 
result  from  assuming  that  the  diffracted  fields  vanish  at  infinity: 

'J  ,U,  *  o(^)  as  ! 


3  U..3U,  *  0(?)  33 
n  d  n  i  s 


(10b) 


The  final  boundary  conditions  are  valid  at  the  upper  and  lower  slit 
Sounder les : 

’V*.zH2<tr-U1(x.z)|z<Q.  for  |  x  |  <a  (Va: 

continuity  of  tangential  H 

'J>(x,2)|2<0,-Ul;x,:)|:<0.  for  jx|<a  ("t! 

€23z,Jl(x,z)|2<0*.cl32'J2;x.z)|z<tj.  for  |x|<a  i'-U. 

continuity  of  tangential 
sl32U2(x,z)|2^0.-e232U1(x,z)|2<>0.  for  |x|<a  (’2b1 

With  the  above  set  of  equations  and  boundary  conditions,  a  unique 
solution  exists  for  the  diffraction  problem. 

Development  of  Integral  Field  Representations 

To  find  the  field,  the  two  dimensional  Green’s  theorem  is  applied,  with 
one  function  given  oy  U(x,o  ar.d  the  other  by  a  standard  Green's  function: 

(7**k*)  u.  ■  -5 ( <-x  * ,  z-z")  : j - 1 ,2,3)  ('3' 

■j  v 

where  ;x,z)  refers  to  a  field  comt  of  interest  and  x’.z'  are  integration 


variables.  Gmce  J 


j  satisfies  I;.  ]).  Green's  theorem  reduces  to: 


♦  IiJ 


0 


'-'orae  and  "aahdacn*:  ' 


0  9) 


Gjtx.zsx'.z’) 


i  ,^12-z'L.l  r 

5177  2a  at! 


y  1 cosf 
m  k 


mir(x»a) 

2a 


Jcos[ 


mir(x**a) 

2a 


]etrn|2. 


where : 

Y  -✓ki-(2I)«  Re(Y  ) >0  Ia(Y  ) >0  (20) 

m  2  2a  a  -  a  — 

With  these  explicit  representations  for  G  and  the  boundary  conditions  on 
both  G  and  U,  the  Integral  representation  for  U  given  in  Eq.  (14)  can  be 
simplified.  Furthermore,  with  the  continuity  conditions  at  the  slit 
boundaries,  U  can  be  written  exclusively  in  terms  of  Ua  and  3ZU„  at  these 
boundaries.  We  find: 

U  (x,z )— /  (— )G,(x,z ; x *  ,b)DUK(x'  )dx*  for  b<z<«  (21a) 

-a  Ei  0 

U ,(x,z)-/  (— )g,(x,z;x'  ,0)DU#(x'  )dx'  for  —  <z<0  *  (21b) 


'J1(x,z)--/a^[GJ(x,z;x'  ,0)0U,(x' )  -U,(x  * )  dz  ,  G  2  (x,z  ;x '  ,  Z’ )  |  z^,  ]dx 

*/a  [Gj( x, z :x, b) ( x ' )-Ub(x ' )92,Gj(x, z ;x' ,z ' ) | z„b-]dx 
for  I x | <a  and  C<z<b 


where  the  boundary  fields  are  defined  by: 


U,  (  X  )  sUj  'v  x  ,Z  )  I  2<q» 

:22a 

DUa(x)53zu1(x,z; | 

:22d 

U(j(x)5Ul(x,z)  | 

(22: 

/ 

X  )  =^2^1  (  X,  Z  i  ,->0- 

'  22  z 

Hence,  the  field  can  te  fzur.d  at  any  zesired  point  once  the  four  unknown 
functions  in  Eqs.  (22)  have  teen  zeteramed.  The  evaluation  of  these 
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functions  is  discussed  below. 

Development  of  Integral  Equations  for  the  3oundary  Functions 
A  set  of  four  integral  equations  is  needed  to  completely  determine  tne 
boundary  functions.  These  equations  result  immediately  from  application  of 
the  continuity  equations  (it)  and  (12)  to  the  field  equations  (21)  at  the  slit 


boundaries.  We  find: 

2ui(x)-UK(x)-/a  (— MG ,  ( x,  b ;x ' , b) DU  (x ' )dx ' 
b  b  _a  e  a  b 

for 

|x|  <a 

(23a) 

Us(x)-/a(^)Gl(x,0:x,,0)D(J,(x')dx' 

for 

|  x  |  <a 

(23b) 

-a  c* 

♦/aO,(x,b;x\a)OUb(x')dx'  for  |x|<a  (230 

jUo(x)-;a[Ga(x,0iX’'.b)Dab(x')-Ub(x1)32lG1(x,0;x',2')|2+B-]dx' 


-/aGj(x,0;x’ ,D)DUa(x’ )dx'  for  |x|<a  (23d. 

-a 

where: 

v 1(x)..lk-(xsla8  •  5C0S9)  |s|<I  ;2»' 

b  2 

These  coupled  equations  for  the  four  boundary  functions  will  be  solved 
numerically  in  the  following  section. 

Approximate  Solution  of  the  Integral  Equations 

To  demonstrate  the  technique  that  will  be  used  to  solve  the  system  in 
(23)*  we  consider  Eq.  (23a).  The  integral  in  this  expression  is  first  divided 
into  N  suo intervals  o*f  width  2a/N,  each  with  a  central  point  xj  given  by: 


Furthermore,  only  the  field  at  the  midpoint  xk  of  each  sub  interval  will  te 
considered.  Hence,  we  have: 

i  S  e,  , 

2U  (x  )-U  (x  )  ■  l  <— )/  G, < x  ,b:x‘,b)DU  (x’)dx' 
b  k  b  k  j«1  e,  CJ  k  b 

J-k 

*(|^)/CkGt(xk,b;x' .b)DUb(x' )dx'  (26) 

The  kth  integral  is  written  independently  since  it  contains  the  point  xk 
corresponding  to  the  singularity  in  G,(xk,b;x' ,b)  as  x’*xk. 

£q.  (26)  can  be  further  reduced  by  applying  the  first  mean  value  theorem 
to  the  N-1  sub  Intervals  Cj  and  the  second  mean  value  theorem  to  the 
subinterval  Ck.  This  yields: 


2VVVV  -j, 


J«X 

‘ft  Uc-^G/.x^.b-.x' ,b)dx' ]ou&(fk)  (27) 

The  above  expression  is  still  exact  for  some  unknown  set  of  N  points 
Zj,  each  within  its  own  corresponding  domain  C*.  However,  an 
approximation  is  now  introduced  to  simplify  the  result,  tn  particular,  when 
the  number  of  subintervals  N  becomes  large,  Zj  can  be  estimated  by  the 
midpoint  of  each  interval: 

*  .  ■  x  ,  £  3  ; 

■J  J 

so  that: 


2WVV 


j  *  2'j.  ,  x  ) 
X,  J  b  j 


:Q  ) 


<-1,2 _ N 


a  *  -  ^  Gt(x,.b;x.,b)  for  >k 

J  NCj  i'  J 


J.k-1 


»»!><  *?t'c,<0‘<Ve  (3°» 

Eq.  (29)  can  be  written  In  matrix  form: 

2uJ-Ub-SIDUb  (3D 

where  the  components  of  the  vectors  are  given  by: 

<Ub5j  "  Ub(ltj)  (32) 

with  similar  expressions  for  Ub,DUb  (and  U0 , DU,  later  on).  The  matrix 
elements  are  found  from  the  expllcity  representation  of  G,  given  in  Eq.  (17a) 

\Zti  -  (^Ix-Jl)  for  j-k  (33a 

3k.k  • 1  ntt  [h«  (ir}  i  ^  i(ir)H#  (ir)n 

(33b 

The  Struve  functions  m0  and  Ht  result  from  an  integral. on  of  the  Hankel 
(1) 

function  H„  about  its  singularity  as  outlined  in  Abramowitz  and  Stegun.7 

Application  of  this  same  technique  to  the  integral  in  Eq.  (23b)  yields  a 
second  matrix  equation: 


s*“  DU. 


where: 


.III  ,a  £,  u(l)  .2k,a,.  , ,  , 


for  *«k 


“  lSff  ("•  (¥>*i  ; 

Note  tr.at  SA  *  - ,  like  S*,  '.as  only  v  independent  elements  and  hence  can  be 


readily  calculated. 


Finally,  the  integrals  in  Eqs.  (23c)  and  (23d)  can  be  similarly 
approximated.  After  considerable  calculation,  we  obtain: 

•  -  *XIDU,*  oIIu0-»IIOUb 
2$,  -  -srIDU,*  °irUb*  ":I3Ub 

wnere: 

II  i  lY»b  i  •  i  j,k  iYmb 

r  -  ■  ■  ■  a  ♦  —  )  —  t  a 

k,j  2HY,  *  N  a£,  Ya  m 


(36a) 

(36b) 

(37a) 


to  solve  tit  is  system  for  tr.e  vectors  U,,  0b,  DU,,  and  DUb  which  represent 

the  values  of  the  previously  unknown  boundary  functions  at  N  distinct  points. 

♦ 

Toward  this  end,  U,  and  Ub  ran  te  eliminated  from  the  system  to  give  a 
single  2Nx2M  matrix  equat::r. ; 


*•  3  »C 


(33) 


where: 
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Since  a  and  C  are  composed  completely  of  known  elements,  the  unknowns  DU9  ar.d 
DUb  can  be  readily  determined.  SuPstitution  into  Eqs.  (36a)  and  (36b)  then 
determines  the  remaining  unknowns  U,  and  'Jjj.  In  other  words,  by  formulating 
the  problem  in  the  language  of  linear  algebra,  a  discrete  solution  has  been 
found  for  the  complete  diffraction  problem.  As  the  number  of  subdivisions 
increases,  tn:s  discrete  solution  should  approach  t.ne  true  continuous 
solution. 

There  are  several  advantages  to  the  method  outlined  nere.  First,  the 
accuracy  of  the  results  can  ce  increased  without  altering  the  basic  algorithm 

i 

oy  increasing  tne  number  of  suo intervals .  Second,  the  validity’  of  the 
solution  is  not  limited  to  tne  narrow  slit  limit.  Finally,  the 
electromagnetic  fields  ar.c  e'e-gy  fl„x  oar.  te  determined  using  the  simple 


numerical  approacn  described  below. 


<3>  •  7-*^  *  E*xw]  «  ^  (Exm#)  -C a i  ; 

FurtherBort,  sine#  the  object  to  be  resolved  will  lie  in  a  plane  parallel  to 
the  screen  containing  the  slit,  we  are  most  interested  in  the  energy  flux 
through  such  a  plane,  and  thus  concern  ouselves  primarily  with  <S2>, 

Lastly,  the  flux  should  be  normalized  to  the  incident  flux  determined  from 
Eq.  (?) .  This  yields: 


sis  ae<Ex"*> 


Determination  of  Transmission  Coefficients 

Although  the  normalized  transmission  coefficient  for  a  particular 
geometry  could  be  calculated  by  integrating  the  above  expression  over  the 
lower  edge  of  the  slit,  an  expression  which  is  more  tractable  numerically  can 
be  obtained  by  manipulating  the  integral  expressions  for  E  and  h.  This  is 
demonstrated  below. 

Integration  of  the  normalized  flux  over  the  slit  tetermines  the 
transmission  coefficient: 


- L- -  i  lim  i  w**i»«ii>cx 

•acose  *  x  x  ‘ 

-a  z-0 


’J3irg  tr.e  exact  integral  forms  of  Ex  and  H  yields: 
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After  several  changes  in 


•  :e-  :f  :  egr  at  id-. ,  these  expressions  reduce 


to: 
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where: 
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Whan  tha  explicit  fora  for  G,  is  inserted,  I,  and  I,  can  da  evaluated 
analytically.  It  is  found  that: 

I,  -  !«in(lc,|x'-x"|> 

I.  •  -|Hi2)(kJ|x'-x”|) 


(47a) 


(476) 


(48a) 

(48b) 


Substitution  of  these  expressions  into  £q.  (46)  and  subsequent  manipulation  of 
DU„  and  DU, *  leads  to  the  final  form  of  the  transmission  coefficient: 


N  N 
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(49: 


where  J„  represents  the  zeroth  order  Bessel  function.  The  elements  (3U#)j 
are  obtained  immediately  when  the  system  is  solved,  and  J0  is  analytic,  so 
that  the  transmission  coefficient  is  easily  obtained. 


NUMERICAL  METHODS 

The  accuracy  and  range  of  /alidity  of  this  model  is  limited  solely  by  the 
number  of  subdivisions  N  -see  in  approximating  the  integral  equations  as  a 
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matrix  equation.  In  parti:-. ar,  1,  in  tne  analog  of  Eq.  (26)  is  roughly 
constant  over  tne  intervals  ..  -_3ric*  the  first  mean  value  theorem  can  oe 

applied  only  if  the  argument  tne  -la.nxel  function  in  1,  is  roughly 


constant.  Hence,  from  Sqs.  ('8),  we  ®uat  have: 

x  >  |i  ;50) 

for  the  approx laat Iona  to  be  valid.  Following  the  example  of  Neerhoff  and 
Mur,  N-30  was  used  to  evaluate  our  Boundary  functions,  so  that  the  results 
could  be  applied  to  find  the  energy  flux  at  distances  greater  than  1/30  of  the 
slit  wldtn  beyond  the  screen.  The  resulting  60x60  matrix  equations  were 
solved  using  Gaussian  elimination.  The  solution  was  checked  for  numerical 
accuracy  by  back  substitution  into  the  matrix  equation,  and  was  found  to  yield 
differences  in  the  components  of  C  of  less  than'one  part  in  10*. 


COMPARISON  OF  TRANSMISSION  COEFFICIENTS 

The  first  objective  of  this  study  was  to  check  the  consistency  of  the  - 
results  by  comparing  the  transmission  coefficients  determined  from  the 
integral  approach  of  Neerhoff  and  Mur  to  those  obtained  by  the  method  of 
Harrington  and  Auckland.  Tcward  this  end,  Eq.  («9)  was  used  to  calculate  T 
for  four  different  slit  widths  and  a  variety  of  scree.-  thicknesses.  The 
results  are  shown  in  Fig.  u. 

From  this  figure,  it  is  seen  that  the  transmission  resonances  of  A / 2 
periodicity  predicted  by  Harrington  and  Auckland  are  reproduced.  Furthermore, 
tne  resonance  positions  are  very  similar,  and  the  peak3  heights  are  in 
agreement  with  the  expression: 

T  •  —  at  resonance  o’: 

K  A 

ieveloped  by  Harrington  aro  lowland  :n  the  narrow  slit  limit.  In  short,  the 
excellent  agreement  between  •o.ese  two  disparate  methods  lends  confidence  in 


RESULTS  AND  IMPLICATIONS 


Near-Field  Scanning  Optical  Microscopy  (NSCM) 

To  determine  the  degree  of  collioation  within  the  near  field,  the 
boundary  functions  corresponding  to  the  initial  resonant  thicknesses  of  Fig.  4 
were  used  to  calculate  the  energy  flux  in  planes  parallel  to  the  screen.  The 
results  are  shown  in  Figures  5a-5d  for  four  different  slit  widths  (A/2.5 
through  A/20).  Froa  these  curves,  the  full  width  at  half  maximum  (FWHM)  of 
the  energy  flux  perpendicular  to  the  screen  was  determined  as  a  function  of 
the  distance  from  the  screen.  These  results  are  in  Figure  6. 

From  the  data  in  Figures  5  and  6,  it  is  possible  to  estimate  the 
feasibility  of  near-field  microscopy  at  optical  wavelengths.  First,  we  note 
that  tne  outward  travelling  radiation  la  collimated  to  exactly  the  slit  width 
at  the  edge  of  the  screen.  This  is  a  direct  consequence  of  the  perfect 
conductivity  assumption,  but  is  likely  to  hold  in  a  realistic  situation 
involving  thin  metal  films  of  sufficient  opacity.  Her.ce,  the  basic  concept  of 
super-resolution  microscopy  remains  valid:  the  light  passing  through  an 
aperture  in  an  opaque  screen  can  be  used  to  define  a  3ub-wavelength  light 
source  which  can  then  be  used  to  generate  an  image  by  scanning. 

Further  analysis  of  Figures  5  and  5  indicates  the  practical  difficulties 
inherent  in  rear-field  microscopy,  due  to  the  divergence  of  the  radiation 
past  the  slit.  As  a  rough  criterion,  the  full  width  at  half  maximum  of  the 
radiation  distribution  is  no  more  than  205  greater  than  the  geometric 
projection  of  the  slit  at  a  ii3tance  of  approximately  .naif  the  slit  width.  In 
the  visible  regime  (\*50CCA;,  t.nis  implies  that  an  aperture  must  be  placed 
within  200A  relative  to  a  surface  in  order  to  obtain  an  image  with  =50CA 
resolution.  If  we  wish  charges  ir  tne  ietected  light  intensity  to  be  due  to 
variations  in  optical  propert.es  cf  tne  sample,  then  tr.e  snarp  decrease  in  the 
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peak  flux  as  a  function  of  distance  places  a  further  constraint  on  the 
precision  of  the  aperture  positioning.  When  scanning  over  low  contrast 
samples  such  as  biological  specimens,  this  suggests  that  a  positioning 
accuracy  of  =20A  must  be  maintained  In  order  to  achieve  a500A  resolution  using 
visible  light. 

These  conclusions  clearly  illustrate  the  techr ical  challenges  inherent  in 
extending  the  microwave  results  of  Ash  and  Hicholls  to  the  optical  regime. 
However,  due  to  recent  developments  in  micropositioning  technology  such  as 
seen  In  the  scanning  tunneling  microscope*  and  the  capacitance  microscope*, 
even  the  strict  requirements  of  near-fleld  optical  microscopy  can  be  met. 
Indeed,  by  employing  microposl cloning  technology,  we  have  already  developed  a 
prototype  microscope,  the  Near-Field  Scanning  Optical  Microscope  (NSOM).1#,l‘ 
With  this  instrument  we  have  obtained  better  than  10G0A  resolution  using  a 
light  source  with  X-5000A,  thus  circumventing  the  conventional  diffraction 
limit.  Fischer1 *  and  Pohl,  at  al.1*  have  also  worked  on  the  problem  of 
near-fleld  optical  imaging,  and  nave  similarly  conduced  that  subwavelength 
resolution  is  achievable.  In  addition  Massey1*  has  investigated  the  infrared 
case,  -where  the  positioning  requirements  are  less  severe  because  of  the  larger 
wavelengths  used. 

Although  half  of  the  slit  width  can  be  used  to  provide  an  estimate  of  the 
dollimation  distance,  it  is  also  clear  from  Figure  6  tr.at,  even  when  the 
distance  is  scaled  to  the  slit  width,  the  radiation  tends  to  diverge  more 
rapidly  as  the  width  decreases  particularly  for  the  smaller  widths.  Hence, 
tr.e  apertures  used  fop  sr.cdd  not  ce  made  so  small  as  to  prohibit 

positioning  of  the  sample  -itrir.  tr.e  collimated  rear-field,  limits  on  signal 
detection  capabilities  also  :e  a  lower  bound  to  tr.e  aperture  size.  In 
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the  other,  an  upper  bound  is  provided  by  the  fact  that  the  aperture  should  oe 
at  least  as  small  as  the  desired  resolution. 

The  existence  of  transmission  resonances  as  a  function  of  screen 
thickness  leads  to  the  question  of  what  effect  the  screen  thickness  has  on  the 
spatial  distribution  of  the  energy  flux.  Presumably,  the  high  transmission  at 
resonance  occurs  when  the  system  efficiently  channels  radiation  from  a  wide 
area  through  the  slit.  In  this  case,  one  might  assume  that  if  the  energy  flew 
is  symmetric  about  the  screen,  then  the  radiation  should  diverge  very  rapidly 
past  the  screen  at  resonance.  Thus,  the  large  signal  strength  associated  with 
resonance  could  only  be  obtained  at  the  expense  of  the  resolution  of  the 
system.  To  test  this  hypothesis,  energy  flux  distribution  curves  similar  to 
Figures  5a-5d  were  developed  for  four  slit  widths  corresponding  to  the 
anti-resonance  positions  in  Figure  u.  By  comparing  the  two  sets  of  data  it 
was  found  that  the  resonance  effect  does  not  influence  the  collimatlon.  To 
verify  this  unexpected  result,  the  FWHM  of  the  energy  flux  was  calculated  for 
the  four  different  slit  widths  as  a  function  of  scree-,  thickness  at  one 
particular  distance  from  t.-.e  screen.  It  was  seen  that,  except  for  very  small 
thicknesses,  the  distr ibut ion  is  indeed  independent  of  thickness.  Thus,  mgr. 
transmission  oar.  oe  achieved  without  a  concurrent  loss  in  -“solution.  In 
-etrospect,  this  result  is  reasonable,  since  t.-.e  symmetry  of  the  prco.em  is 
disrupted  Oy  the  presence  of  t.-.e  initial  and  reflected  fields  in  addition  to 
tne  diffracted  field  on  or.e  3ide  (£q.  9).  i3  the  thiowr.ess  changes,  the 

♦  j  *  ♦ 

fields  ja  and  'J,  change  only  .-.  magnitude,  out  tne  field  J l  changes  m 
distribution  as  wen  eir.ee  .-.valves  t.-.e  sum  of  - 3  with  tr.e  unchanging 
f  .e  Ids  ^  ar.d  Jr .  it  -esd-.ar :  e .  tne  distribution  of  J .  leads  to  a 
channeling  of  t.-.e  radiatio- ,  mo  o-e  disomouoion  of  * , 


remains  unaffected. 


3-bwaveler.gtn  Photolithography 


The  near-field  colliaation  phenomenon  can  also  be  exploited  to  produce 
subwaveleng tit  structures  via  pnotolithograpny  in  the  visible  or  near  UV 
regime.  However,  because  the  collimation  range  is  so  short  at  optical 
frequencies,  the  pattern  in  a  mask  can  be  reproduced  with  high  fidelity  in 
photoresist  only  if  the  mask  is  in  close  proximity  to  the  resist  and  only  if 
the  resist  Itself  is  very  thin.  Again  referring  to  Figure  6,  the  mask  should 
be  roughly  2Q0A  or  less  away  from  the  furthest  side  of  the  resist  in  order  to 
produce  500A  structures.  This  may  be  achievable  by  suitably  thinning 
conventional  resists,  but  should  certainly  be  possible  using  recently 
developed  Langmuir-Blodgett  film  technology.1* 

This  method  of  super-resolution  pattern  generation  has  already  been 
demonstrated  oy  Fischer  and  Zingsheim. l*  Using  a  conventional  resist  and  a 
blue  light  source,  they  obtained  100QA  resolution,  and  with  a  photobleachir.g 
dye  layer  they  obtained  50C-’00A  resolution.  We  are  currently  working  to 
extend  these  results  using  sharper  masks  formed  by  electron  beam  lithography. 

The  model  studied  here  is  more  applicable  to  lithography  than  to 
microscopy,  since  the  masks  used  in  the  production  of  integrated  circuits  are 
generally  formed  from  multiple  slit-like  structures,  -hereas  the  apertures 
used  for  two  dimensional  'JaCM  must  be  approximately  round.  However,  the  mere, 
is  restrictive  in  that  the  magnetic  field  .3  assumed  to  oe  aiir.ged  along  tr.e 
slit.  This  leads  to  the  existence  of  propagating  modes  within  tr.e  slit  'and 
correspondingly  large  transmission  coefficients),  out  for  otr.er  polar izatiors 
the  modes  may  be  evanescent.  Trus,  if  polarized  light  is  used  for 

t 

subwavelength  pnotolithogr aory ,  tr.e  exposure  of  the  resist  may  depend  on  tr.e 
orientation  of  tr.e  mask  as  -e..  is  or.  tr.e  sizes  of  tr.e  features  incorporated 
therein.  To  avoid  thi3  .-00. anted  :  1  l-mmat  ion  should  oe  ^3ed. 
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The  region  Between  the  mask  and  the  substrate  underlying  the  realat  can 
alao  form  a  two  dlmenaional  waveguide  'which  aupporta  propagating  modes.1*  If 
monochromatic  illumination  ia  uaed,  theae  modes  can  produce  a  aeries  of  A/2 
spaced  fringes  in  the  developed  resist  t*Uch  are  analagous  to  the  transmission 
resonances  shown  in  Figure  u .  However,  if  polychromatic  illumination  is  used, 
the  field  maxima  under  the  mask  should  occur  at  different  points,  so  that  the 
fringes  are  eliminated  By  destructive  interference.  For  example,  a  Xe  arc 
lamp  would  be  preferrable  to  a  laser  or  even  to  a  Hg  arc  lamp  because  it  has  a 
smooth,  broadband  output  in  the  visible  and  near  UV. 

High  Density  Optical  Information  Storage 

The  NSOM  concept  can  be  combined  with  the  principle  of  subwavelength 
photolithography  to  produce  high  density  read-write  optical  memories.  To 
accomplish  this,  the  photoresist  must  be  replaced  by  a  bistable  optical 
material;  that  is,  a  material  whose  optical  properties  (e.g.,  absorption)  can 
flip  reversibly  between  two  states  based  on  the  intensity  and/or  the  frequency 
of  the  radiation  employed.  The  light  passed  through  an  aperture  is  then  used 
to  change  a  small  section  of  the  surface  to  either  t re  "on”  or  "off"  state. 

The  information  in  many  such  sections  can  then  be  read  by  using  the  same 
aperture  with  a  different  incident  intensity  and  frequency  as  an  NSOM  prcce. 
The  bit  3ize  is  limited  only  by  the  size  of  the  aperture.  Assuming  that  t re 
bits  are  5CCA  in  diameter  and  are  separated  by  50CA,  'l13  bits  (more  than  ' 
gigabyte)  of  information  can  be  stored  in  1  omJ. 

Near-Field  Scanning  Acoustic  Microscopy  (N5AM ) 

The  above  analysis  of  tie  rear-field  is  also  directly  applicable  to  N5AM, 
Near-Field  Scanning  .Aeoust.:  vi croscopy .  In  the  acoustic  case,  the  magnetic 
field  'J  is  replaced  by  the  cressure  3 .  This  field  still  satisfies  the 
H.eimnoltz  equation  if  tne  ex:.tatior.  .3  rarmcmc: 
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where  0  is  the  density  of  the  fluid  (analagous  to  the  index  of  refraction). 
The  velocity  field  u  is  derivable  from  P  just  as  the  electric  field  E  is  from 
the  magnetic  field,  ,  except  that  it  is  longitudinal  to  the  direction  of 
travel  rather  than  transverse  to  it: 

3  *  £p  (54) 

Furthermore,  if  the  assumption  of  perfect  conductivity  is  replaced  by  the 

assumption  that  the  screen  is  perfectly  rigid,  then: 

3  P-0  at  the  screen  (55) 

n 

Comparing  Eqs.  (52)  and  (55)  to  Eqs.  (3)  and  (9)  it  is  seen  that  P  satisfies 

exactly  the  same  differential  equation  and  boundary  conditions  as  does  U. 

Hence,  the  above  numerical  solution  for  the  magnetic  field  immediately  yields 

the  pressure  in  the  acoustic  case,  and  the  results  concerning  the  collimatlon 

distance  and  the  feasibility  of  super-resolution  imaging  remain  valid. 

NSAM  has  the  potential  to  extend  the  resolution  of  conventional, 

diffraction  limited  acoustic  microscopy* *  in  the  same  manner  that  NSOM  can 

extend  optical  microscopy.  However,  NSAM  has  two  additional  advantages. 

First,  acoustic  microscopy  i3  limited  in  resolution  by  attenuation  within  tne 

medium  supporting  the  wave.  Since  one  could  obtain  \/’0  resolution  using  the 

near-field  technique,  the  frequency  of  tne  incident  acoustic  wave  can  be 

reduced  by  an  order  of  magnitude  with  no  concurrent  loss  of  resolution.  The 

/ 

d3  attenuation,  proport icr.al  to  the  frequency  squared,  can  thus  be  reduced  by 
a  factor  of  'TO  without  loss  :f  -esolucion  if  an  NSAM  system  with  a  A / i 0 


aperture  is  employed.  Second,  a  cylindrical  acoustic  waveguide,  such  as  an 
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aperture  in  a  thick  rigid  screen,  supports  at  least  one  propagating  mode  no 
matter  how  narrow  the  constriction.  Hence,  attenuation  of  the  signal  within 
the  aperture  should  be  small  compared  to  the  optical  case. 

Topographical  Mapping  and  Mask  Alignment 

The  rapid  decrease  in  the  transmitted  intensity  past  a  slit  can  also  be 
used  to  advantage  in  two  ways.  First,  a  nondestructive  topographical 
instrument  could  be  constructed  by  collecting  light  passed  through  a  slit  or 
aperture  and  reflected  back  through  the  aperture  by  an  object  within  the 
near-field.  As  the  aperture  is  scanned  over  an  optically  uniform  surface,  the 
intensity  of  the  reflected  light  would  represent  changes  in  the  aperture  to 
object  separation  due  to  variations  in  the  height  of  the  object.  Second,  this 
same  reflection  scheme  with  a  feedback  mechanism  could  be  used  to  position  a 
lithographic  mask  with  respect  to  a  substrate. 

The  analysis  of  the  near-field  diffraction  by  a  slit  could  be  applied  to 
a  number  of  other  problems  as  well  ranging  from  electromagnetic  Interference 
between  close  packed  semiconductor  devices  to  the  shielding  of  an  environment 
from  RF  noise. 

CAVEATS 

The  Limitations  of  the  above  analysis  must  oa  considered  before  the 
results  are  used  to  apply  near-field  phenomena  to  experimental  situations. 
First,  it  3hould  oe  noted  tnat  the  slit  model  was  chosen  more  for  its 
mathematical  tractadllity  than  its  correspondence  to  reality.  For  example, 
the  metal  films  which  can  be  -sed  to  form  a  thin  screen  have  finite 
conductivities  and  ne/sce  are  net  completely  opaque.  This  would  increase  the 
noise  in  an  VSCM  system  Cue  tc  tr.e  passage  of  a  small  fraction  of  the  light  :r. 
the  region  away  from  the  aperture.  Furthermore,  the  effective  diameter  of  an 


aperture  would  increase,  since  energy  could  propagate  below  the  surface  of  the 
metal  wall  forming  the  aperture. 

The  model  considered  above  Is  also  limited  to  the  situation  where  the 
incident  magnetic  field  is  aligned  along  the  slit.  For  other  polarizations, 
the  transmitted  intensity  may  be  radically  different,  but  the  degree  of 
coll imat ion  should  be  roughly  the  same. 

In  interpreting  the  results,  it  should  be  remembered  that  the  full  width 
of  half  maximum  of  the  outgoing  flux  provides  only  a  rough  estimate  of  the 
resolution  to  be  expected  in  a  NSOM  system.  Of  additional  importance  is  the 
sharpness  of  the  distribution  below  the  slit,  as  is  evident  in  Figures  5a-5d. 
Using  this  additional  criterion,  it  is  seen  that  the  FWHM  probably  overstates 
the  resolution,  and  that  it  may  be  necessary  to  place  the  aperture  closer  tnan 
half  the  width  in  order  for  the  resolution  to  approach  the  aperture  size. 

Another  potential  problem  concerns  the  effect  of  the  object  (for  NSOM)  or 
the  resist  and  substrate  (for  lithography)  on  the  distribution  of  radiation 
emanating  from  an  aperture.  Clearly,  any  such  objects  placed  within  the 
near-field  change  the  boundary  conditions  and  hence  the  solution  of  the  entire 
problem.  The  most  significant  changes  should  occur  for  the  case  of  conducting 
surfaces,  but  in  all  situations  one  might  assume  that  the  substrate  should 
cause  the  radiation  to  spread  more  rapidly,  since  the  screen  and  the  surface 
can  act  as  a  waveguide  to  funnel  radiation  away  from  the  aperture. 

Although  the  principle  of  r.ear-field  imaging  can  be  investigated  using 
tne  infinite  slit  model  stymied  .mere,  a  circular  aperture  must  be  used  for 
dimensional  imaging  in  *liiv  applications.  In  terms  of  the  transmitted 
intensity,  this  can  nave  s.gr.fioant  consequences:  tne  slit  m  a  thick  screen 
acts  as  a  waveguide  whic.r  ran  support  a  propagating  ncce,  wnereas  a  round 
aperture  in  a  thick  screen  ran  a.cccrt  crly  evanescent  nodes  if  tne  aperture 


is  considerably  smaller  than  a  wavelength,  r urthermore,  there  is  an  additonal 
reduction  in  transmission  due  to  reflections  at  the  screen  which  occur  even  in 
the  simpler  case  of  an  Infinitely  thin  screen. 

The  geometry  of  an  aperture  closed  in  two  dimensions  has  been  considered 
by  many  researchers.  Massey1*  studied  a  square  aperture  in  a  thin  screen 
using  Fourier  optics  and  considered  the  feasibility  of  super-resolution 
microscopy  and  lithography.  The  qualitative  conclusions  of  the  work 
concerning  the  extent  of  collimatlon  are  consistent  with  the  more  extensive 
distribution  graphs  in  Figures  ua-4d,  although  the  approach  used  may  not  be 
valid  in  the  near-fleld.  More  recently,  Leviatan1*  has  determined  the 
distribution  of  the  magnetic  field  transmitted  by  a  circular  aperture  in  a 
thin  screen  using  the  method  of  moments.  The  collimatlon  results  were  also' 
found  to  be  similar  to  the  slit  geometry  studied  here.  The  more  realistic 
case  of  a  round  aperture  in  a  thick  screen  has  been  considered  by  McDonald21, 
although  only  the  degree  of  coupling  between  the  two  sides  of  the  screen  was 
investigated  rather  than  the  distribution  of  radiation  in  the  near-field. 

The  results  of  McDonald  can  be  used  to  obtain  a  -ough  estimate  of  the 
minimum  aperture  size  which  can  be  used  for  practical  NSOM  imaging.  The  least 
attenuated  mode  within  a  round  aperture  was  found  to  :e  the  TE x  t  mode,  for 
wnich  the  energy  decays  at  a  rate  of: 

_  -2 *1.841  d/a 

-  -  i,e  ,oo  , 

where  d  is  the  screen  thickness  and  a  is  the  aperture  radius.  However,  the 
energy  decay  in  the  screen  it  a  finite  rate  as  well,  «r. icr.  is  given  by: 

-  -3/s 

-  *  -3 9  '  ' 

where  5  is  the  extinction  .e-gtr.  ;r.  the  screen.  The  metal  with  the  .argest 


opacity  in  the  visible  :s  a .  .  - .  - for  -r..;r,  5 -65  A  wr.en  \*50C1A.22  When  the 
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attenuation  due  to  the  waveguide  effect  exceeds  the  attenutation  in  the  metal, 
the  contrast  between  the  aperture  and  the  surrounding  screen  becomes 
insufficient  for  super-resolution  applications.  Hence,  the  minimum  aperture 
diameter  is  given  by: 

2a  -  4*1.841.6  (58) 

For  At,  this  yields  a  diameter  of  *480A,  which  explains  our  choice  of  a  goal 
of  500A  for  NSOM  imaging.  It  should  be  noted  that  this  is  a  very  approximate 
result,  since  the  finite  conductivity  of  the  screen  should  decrease  the 
waveguide  attenuation  below  the  value  given  in  Eq.  (56)  due  to  the  larger 
effective  aperture  size  as  discussed  above. 

In  a  strict  sense,  the  resolution  of  an  NSOM  system  should  indeed  by 
limited  to  the  aperture  size,  since  two  or  more  small  structures  will  not  be 
completely  resolved  if  their  separation  is  less  than  the  aperture  diameter. 
However,  it  should  be  possible  to  determine  the  position  of  the  edge  of  a 
structure  larger  than  the  aperture  (i.e.,  determine  its  position)  to 
considerably  less  than  the  aperture  size.  To  accomplish  this,  the  intensity 
information  obtained  when  the  aperture  is  partially  occluded  by  the  edge  would 
be  deconvoluted  to  find  the  edge  position. 

Finally,  although  it  is  easiest  to  understand  near-field  microscopy  in 
the  senario  where  the  aperture  is  used  for  form  a  subwavelength  sized  light 
source,  thi3  is  by  no  means  the  only  possible  mode  for  (ISCM  viewing.  In 
addition,  the  light  transmitted  through  the  sample  cculd  be  collected  2Z  :ne 
aperture  in  the  near-field.  In  this  arrangement,  only  the  light  from  the 
region  directly  below  /the  aperture  is  detected.  Although  the  slit  model 
suggests  that  thi3  geometry  may  yield  a  lower  resolution  for  a  given  slit  to 
object  separation  due  to  t re  -atiatio.n  scooping  effect  at  resonance,  the 
resolution  3hould  not  be  affected  .n  tre  :ase  of  a  circular  aperture  since  all 
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.•nodes  are  evanescent  and  no  resonances  occur.  Finally,  the  light  could  be 
both  transmitted  by  and  collected  with  the  aperture,  thus  forming  a  reflection 
NSOM  instrument.  This  system  has  two  principle  advantages:  even  thick, 
opaque  samples  could  be  scanned,  and  the  resolution  and  contrast  may  be  higher 
than  in  the  transmission  modes,  since  the  neer-field  colllmation  phenomenon  is 
used  to  advantage  not  once  but  twice. 

CONCLUSIONS 

The  energy  flux  transmitted  by  an  Infinite  slit  in  a  thick  conducting 
screen  has  been  determined  using  the  method  of  Neerhoff  and  Mur.  The  model  is 
limited  only  in  the  number  of  subdivisions  taken  to  solve  the  resulting 
integral  equations  by  numerical  methods.  Transmission  coefficients  were  found 
to  be  Identical  to  those  calculated  by  the  alternative  method  of  Harrington  ‘ 
and  Auckland.  The  flux  perpendicular  to  the  screen  is  indeed  collimated  to 
the  slit  width  within  the  very  near-field,  but  diverges  sufficiently  rapidly 
so  that  an  object  must  be  placed  closer  than  half  the  3lit  width  in  order  to 
generate  an  image  with  a  resolution  approaching  this  width.  In  addition  to 
near-fleld  microscopy,  other  applications  such  as  subwavelength  photolitho¬ 
graphy  and  optical  information  storage  were  considered.  Many  physical 
phenomena  which  nay  be  associated  with  the  optical  near-field  were  discussed 
qualitatively  in  conjunction  with  these  applications.  Finally,  the 
Limitatior.3  of  the  slit  model  were  viewed  in  the  context  of  its  application  to 
the  design  and  construction  of  a  practical  Near-Field  Scanning  Optical 
Microscope. 
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figure  captions 


Figure  1.  Geometry  of  the  Near-Field  Scanning  Optical  Microscope. 

Figure  2.  Model  for  diffraction  by  a  slit  in  a  thick  screen. 

Figure  3.  Regions  of  integration  used  for  the  boundary  value  problem. 

Figure  «.  Transmission  coefficients  for  1-5QQ0A  radiation  as  a  function  of 
screen  thickness  for  slits  of  1/2.5,  1/5,  1/10,  and  1/20  widths. 

Figure  5.  Distributions  of  S2,  the  component  of  the  energy  flux  normal  to 
the  screen,  transmitted  by  silts  of  width:  a)  1/20;  b)  1/10;  c) 
1/5;  d)  1/2.5  for  1-5000A  radiation. 

Figure  6.  Rate  of  spreading  of  1-5000A  radiation  past  the  screen  as 

determined  by  the  full  width  at  half  maximum  of  Sz  for  slits  of 

1/2.5,  1/5,  l/’O,  and  1/20  widths. 
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ABSTRACT 


A  new  method  for  high  resolution  imaging,  near  field  scanning  optical 
microscopy  (NSOM),  has  been  developed.  The  concepts  governing  this  method  are 
discussed,  and  the  technical  challenges  encountered  in  constructing  a  working 
NSOM  instrument  are  described.  Two  distinct  methods  are  presented  for  the 
fabrication  of  well  characterized,  highly  reproducible,  sub-wavelength 
apertures.  A  sample  one  dimensional  scan  is  provided  and  compared  to  the 
scanning  electron  micrograph  of  a  test  pattern.  From  this  comparison,  a 
resolution  of  better  than  1500A  (i.e.,  — \ / 3 . 6 )  is  determined,  which  represents 
a  significant  step  towards  our  eventual  goal  of  500A  resolution.  Fluorescence 
has  been  observed  through  apertures  smaller  than  600A  and  signal-to-noise 
calculations  demonstrate  that  fluorescent  imaging  should  be  feasible.  The 
application  of  such  imaging  is  then  discussed  in  reference  to  specific 
biological  problems.  The  NSOM  method  employs  non-ionizing  visible  radiation 
and  can  be  used  in  air  or  aqueous  environments  for  non-destructive 
visualization  of  functioning  biological  systems  with  a  resolution  comparable 
to  that  of  scanning  electron  microscopy. 
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INTRODUCTION  ! 

( 

Resolving  components  spatially  in  biological  systems  presents  a  classical  > 

* 

biophysical  problem.  Most  imaging  techniques  presently  in  use  are  | 

fundamentally  limited  by  the  wavelength  of  the  radiation.  In  this  paper  we  I 

I 

present  our  development  of  a  unique  new  imaging  method  that  is  independent  of 
the  radiation  wavelength.  This  method  permits  the  non-destructive  imaging  of 
surfaces  within  their  native  environment,  with  a  resolution  comparable  to 
scanning  electron  microscopy.  We  have  called  this  advance  in  biophysical 
technology,  near  field  scanning  optical  microscopy  (NSOM) . 

The  resolution  limits  for  this  new  technique  fall  between  the  extremes 
defined  by  light  microscopy  (2500A)  and  fluorescence  energy  transfer 
techniques  (  =  40A  to  -80A).  Using  visible  wavelengths,  a  conservative  estimate 
of  the  resolution  is  approximately  500A .  Thus,  it  is  likely  that  this 
technique  will  make  an  important  contribution  to  the  investigation  of 
macromolecular  assemblies  such  as  lipid  microdoraains ,  receptor  clusters,  and 
photosynthetic  reaction  center  aggregates.  Furthermore,  unlike  other 
techniques  at  this  resolution,  near  field  imaging  will  be  able  to  follow  the 
temporal  evolution  of  these  macromolecular  assemblies  in  living  cells. 

Therefore,  NSOM  has  the  potential  to  provide  both  kinetic  information  and  high 
spatial  resolution. 

The  Conceptual  Basis  of  Near  Field  Scanning  Optical  Microscopy  (NSOM) 

The  fundamental  principle  underlying  the  NSOM  concept  is  outlined  in 

Figure  1,  where  visible  light  is  depicted  as  being  normally  incident  upon  a 

conducting  screen  containing  a  small  (sub-wavelength)  aperture.  Because  the 

screen  is  completely  opaque,  the  radiation  emanating  through  the  aperture  and 

into  the  region  beyond  the  screen  is  first  collimated  to  the  aperture  size 

rather  than  to  the  wavelength  of  the  radiation  employed.  This  occurs  in  the 
near  field  regime.  Eventually  the  effect  of  diffraction  is  evidenced  as  a 


marked  divergence  in  Che  radiation,  resulting  in  a  pattern  that  no  longer 
reproduces  the  geometrical  image  of  the  aperture.  This  occurs  in  the  far 
field  regime. 


To  apply  the  colliraation  phenomenon,  an  object,  such  as  a  cell  membrane, 
is  placed  within  the  near  field  region  relative  to  an  aperture.  In  this  case, 
the  aperture  acts  as  a  light  source  whose  size  is  not  limited  by  the 
considerations  of  geometrical  optics.  The  light  source  can  be  scanned  over 
the  object,  and  the  detected  light  can  be  used  to  generate  a  high  resolution 
image.  Because  the  resolution  is  dependent  upon  the  aperture  size  rather  than 
the  wavelength,  it  should  *be  possible  to  obtain  500A  or  better  resolution 
if  a  sufficiently  small  aperture  is  used.  These  are  the  essential  features  of 
the  NSOM  technique,  which  can  be  applied  in  air  or  aqueous  environments  using 
non-ionizing,  visible  radiation.  Thus,  raacroraolecular  assemblies  in 
functioning  biological  systems  can  be  imaged  at  high  resolution. 

Comparison  With  Other  Optical  Methods 

It  is  important  to  distinguish  NSOM  from  more  conventional  diffraction 
limited  techniques.  For  example,  although  images  have  been  obtained  of 
isolated  250A  structures  using  video  enhanced  contrast  optical  microscopy,1 
the  images  of  these  structures  are,  in  the  best  of  cases,  limited  by 
diffraction  to  sizes  of  1000-2000A .  Hence,  with  this  method  it  is  difficult 
to  determine  the  true  sizes  of  the  observed  features.  Furthermore,  two  or 
more  structures  will  not  be  individually  resolved  if  their  respective  Airy 
disks  overlap  (as  determined  by  the  Rayleigh  criterion).  In  comparison,  the 

t 

NSOM  method  should  generate  images  which  indicate  the  true  sizes  of  assemblies 
even  when  these  assemblies  are  separated  by  much  less  than  the  wavelength. 
Historical  Background 

As  far  back  as  a  decade  ago,  the  principle  of  super-resolution  microscopy 
was  demonstrated  at  microwave  frequencies  (\*3cm)  by  Ash  and  Nicholis.  In 
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their  pioneering  experiment,  a  grating  of  0.5mm  periodicity  was  imaged  with  an 

effective  resolution  of  X/60.  However,  until  we  reported  our  initial  results 
3 

in  1983  ,  there  were  no  published  attempts  to  extend  this  technique  to  the 

visible  region  of  the  spectrum.  This  is  quite  understandable,  because  the 

minute  physical  dimensions  of  the  optical  near  field  demand  aperture 

fabrication  and  micropositioning  technologies  on  a  nanometer  scale. 

Furthermore,  it  is  not  immediately  obvious  that  the  results  of  the  microwave 

experiment  could  be  extended  to  tne  visible  regime,  because  the  finite 

conductivity  of  thin  metal  films  has  to  be  considered  in  the  design  of  screens 

for  submicron  apertures,  whereas  the  comparatively  thicker  metal  screens  used 

in  the  microwave  case  are  completely  conducting  and  thus  opaque. 

Since  our  original  report,  several  additional  publications  on  near  field 

4  11  8 

imaging  technology  have  appeared.  Fischer  has  produced  an  interesting 

series  of  results  by  scanning  a  sub-wavelength  aperture  over  a  second,  larger 
aperture,  but  these  results  are  difficult  to  interpret  for  a  number  of 
reasons.  For  example,  the  opacity  of  the  metal  films  employed  was  not  large, 
so  that  the  apertures  were  poorly  defined  and  much  stray  light  was  transmitted 
through  the  aperture  screens.  In  addition,  coherent,  monochromatic 
illumination  was  used  at  grazing  incidence,  so  that  a  series  of  standing  waves 
may  have  been  generated  when  the  polarization  of  the  electric  field  was 
perpendicular  to  these  screens. 

A  near  field  imaging  system  for  use  in  the  far  infrared  is  also  currently 
being  developed  by  Massey.*0  Although  this  sytera  may  find  many  applications 

t 

(such  as  in  the  detection  of  heat  transport  on  a  microscopic  scale),  it  will 
not  approach  the  resolution  capabilities  of  NSOM  and  can  best  be  viewed  as  a 
complementary  imaging  technique. 

Finally,  Pohl  et  al.**  have  developed  a  system  for  super-resolution 
microscopy  that  shares  some  of  the  same  features  as  our  NSOM  instrument.  They 


cal L  their  technique  opticaL  stethoscopy.  This  is  somewhat  of  a  misnomer, 
as  the  analogy  to  a  medical  stethoscope  is  not  strictly  valid.  This  can  be 
understood  by  treating  the  optical  aperture  and  the  acoustic  stethoscope  as 
classical  waveguides.  In  the  optical  case,  a  consideration  of  the  appropriate 
boundary  conditions  indicates  that  all  modes  are  evanescent  (decaying)  within 
the  aperture.  Indeed,  only  the  finite  thickness  of  the  aperture  screen 
permits  the  flow  of  any  energy  at  all.  A  different  set  of  boundary  conditions 
must  be  invoked  in  the  acoustic  case.  These  conditions  lead  to  the  existence 
of  a  single  propagating  plane  wave1 3  and  explain  why  a  20Hz  (\=17ra)  sonic 
disturbance  can  be  detected  through  a  long  stethoscope  tube  only  a  centimeter 
in  diameter. 

At  the  experimental  level,  the  sizes  and  the  structure  of  the  apertures 
used  by  Pohl  et  al.11  were  not  characterized.  In  addition,  the  nature  of 
their  manufacture  presents  considerable  challenges  in  the  attainment  of  repro¬ 
ducibility.  To  demonstrate  250A  resolution  using  optical  techniques,  not  only 
must  the  apertures  be  well-characterized,  but  it  is  also  essential  to  have 
well-characterized  test  structures.  Neither  of  these  criterion  was  met  in 
Pohl  et  at . ' s  paper. 

In  short,  initial  attempts  to  implement  the  near  field  scanning  concept 

attest  to  the  difficulty  of  the  NSOM  technique.  To  further  appreciate  the 

technical  challenges  inherent  in  this  form  of  microscopy,  we  have  performed 

model  calculations  to  obtain  a  first  approximation  to  the  pattern  of  radiation 

in  the  near  field.  These  calculations  (to  be  published)  investigate  the 

/ 

transmission  of  light  through  a  slit  of  infinite  length  in  a  screen  of  finite 
thickness,  and  were  based  on  the  powerful  Green's  function  technique  of 
Neerhoff  and  Mur.  In  particular,  the  component  of  the  time  averaged 
Poynting  vector  perpendicular  to  the  screen  was  calculated  in  the  region  just 
beyond  the  screen.  The  results  are  shown  in  Figure  2,  where  the  divergence  of 
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the  energy  flux  is  pLotted  as  a  function  of  the  distance  from  the  slit.  These 
results  demonstrate  that  the  radiation  remains  collimated  to  a  distance  of  at 
least  half  the  slit  width.  Furthermore,  they  also  indicate  that  the  extent  of 
the  near  field  increases  with  the  slit  width.  Finally,  these  near  field 
energy  flux  calculations  exhibit  a  close  to  exponential  decrease  in  intensity 
with  increasing  distance  from  the  screen.  This  suggests  that  rigid  stability 
requirements  will  need  to  be  maintained  in  the  z  direction,  as  will  be 
discussed  more  fully  in  the  following  section.  Of  course,  significant 
differences  can  be  expected  between  the  near  field  radiation  pattern  of  a  slit 
and  that  of  a  circular  aperture,  but  the  simpler  slit  geometry  should  give  an 
indication  of  the  pattern  to  be  expected  in  the  aperture  case. 


TECHNICAL  CHALLENGES  TO  BE  ADDRESSED 

To  construct  a  super-resolution  microscope,  several  practical 
considerations  must  first  be  addressed.  These  considerations  fall  into  two 
broad  categories:  raicropositioning  and  enhancement  of  the  signal  to  noise 
ratio. 


Microposit ionint 


The  theoretical  analysis  mentioned  above  indicates  the  degree  of  control 


required  over  the  aperture  to  object  separation.  For  example,  the  radiation 
emanating  from  a  500A  wide  slit  remains  collimated  to  an  approximate  distance 
of  250A.  Furthermore,  because  the  intensity  in  the  near  field  is  so  strongly 


dependent  upon  distance,  the  separation  must  be  maintained  with  a  precision  of 
about  20A .  To  satisfy  such  strict  positioning  requirements,  several  potential 


problems  must  be  successfully  confronted.  These  are  discussed  below. 


Vibration  Isolation 

In  a  typical  laboratory  environment,  building  vibrations  exist  at  fre¬ 
quencies  as  low  as  2  to  4  Hz,  although  most  of  the  vibrational  energy  is  in 
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the  5  to  30  Hz  range.15  In  our  environment,  we  have  measured  displacements  of 
up  to  111  in  this  frequency  range.  If  this  energy  were  to  be  transferred  with¬ 
out  attenuation  to  the  critical  components  of  an  NSOM  system,  and  if  the  aper¬ 
ture  and  object  were  not  held  together  with  sufficient  rigidity,  then  their 
relative  displacement  could  vary  by  2  to  3  orders  of  magnitude  greater  than 
the  precision  required.  To  preclude  this  possibility,  it  is  necessary  to  de¬ 
sign  an  isolation  system  which  attenuates  both  horizontal  and  vertical  vibra¬ 
tions  from  the  floor  at  frequencies  above  a  few  Hertz.  In  addition,  such  a 
system  must  employ  a  damping  mechanism  to  insure  that  any  vibrations  which  are 
transmitted  to  the  microscope  are  transient  in  nature.  This  damping  mechanism 
is  also  needed  to  reduce  the  vibrations  induced  by  the  motion  of  the  sample 
stage  itself  in  the  course  of  a  scan.  Finally,  an  effective  isolation  system 
must  also  shield  the  NSOM  instrument  from  acoustic  vibrations,  which  fall  in 
the  20  Hz  to  20  kHz  range. 

Thermal  Drift 

Although  the  thermal  expansion  coefficients  of  the  various  components  of 
an  NSOM  system  will  vary  widely  with  composition  and  size,  it  can  be  expected 
that  these  components  will  expand/contract  roughly  0.1  to  1.0  microns  for 
every  1 *C  increase/decrease  in  temperature.  Hence,  the  differential  rate  of 
expansion  of  the  aperture  relative  to  an  object  could  present  a  serious 
obstacle  cowards  obtaining  the  required  raicropos it ioning  capabilities  for 
super-resolution  microscopy.  This  obstacle  can  be  surmounted  through  two 
lines  of  attack:  either  the  instrument  can  be  designed  so  as  to  insure  that 
the  thermal  expansion  of  the  aperture  relative  to  the  object  is  rather  small, 
or  else  stringent  control  can  be  maintained  over  the  temperature  of  the  entire 
apparatus.  Ideally,  both  of  these  methods  should  be  used. 
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Aperture  and  Sample  Translation 

As  a  final  raicropositionmg  requirement,  it  is  necessary  to  translate  the 
aperture  and  sample  accurately  relative  to  one  another  in  all  three  spatial 
dimensions.  As  mentioned  above,  the  translating  system  should  position  the 
aperture  with  a  precision  of  roughly  20A  in  the  z  direction  (perpendicular  to 
the  sample  plane).  However,  it  is  equally  necessary  to  obtain  accurate 
positioning  in  the  x  and  y  directions,  as  one  factor  which  limits  the 
resolution  of  an  NSOM  system  is  the  size  of  the  steps  taken  while  scanning. 

In  particular,  the  step  size  should  be  shorter  than  half  of  the  desired 
resolution. 16  Hence,  a  system  designed  for  500A  resolution  must  include 
positioning  control  of  better  than  250A  in  the  x  and  y  directions. 

Signal  Considerations 

The  light  transmitted  through  a  subraicron  aperture  is  weak  enough  to 
require  the  use  of  sensitive  detection  electronics.  A  good  detection  system 
becomes  increasingly  important  when  it  is  considered  that,  in  the  interest  of 
producing  high  resolution  images  rapidly,  the  aperture  must  scan  quickly  over 
the  object,  so  that  the  time  available  at  each  point  to  accumulate  the 
requisite  counting  statistics  is  rather  brief.  Weak  signals  can  particularly 
be  expected  in  biological  applications,  where  contrast  differences  in  a 
specimen  can  be  quite  small.  In  such  applications,  a  variety  of 
non-destructive  methods  would  be  needed  to  increase  contrast,  ranging  from  the 
computer  enhancement  of  images  to  the  fluorescent  labeling  of  specimens. 

t 

Noise  Considerations 

Since  rather  weak  signals  are  expected,  one  obvious  source  of  noise  is 
that  of  statistical  noise.  The  significance  of  this  source  can  be  reduced  by 
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increasing  the  period  of  data  collection,  but  this  is  done  at  the  expense  of 
scan  speed.  Systematic  noise  can  also  be  expected  from  the  detection  system 
in  the  form  of  dark  noise  and  noise  from  the  preamplifier  and  discriminator 
electronics.  Noise  will  also  be  introduced  if  there  is  any  variation  in  the 
light  intensity  used  to  illuminate  the  sample  or  aperture.  Furthermore,  since 
the  aperture  screen  is  not  completely  opaque,  any  light  which  is  transmitted 
through  this  screen  (rather  than  through  the  aperture)  will  contribute  to  the 
systematic  noise. 

There  are  also  sources  of  noise  related  to  the  issue  of 
raicropositioning.  First,  unwanted  variations  in  the  xy  positioning  of  the 
sample  relative  to  the  aperture  can  either  cause  an  apparent  distortion  in  the 
image  (for  variations  long  compared  to  the  scan  time)  or  else  can  smear  the 
contrast  and  resolution  information  at  a  point  (for  variations  short  compared 
to  the  data  collection  time  af.  that  point).  Of  even  more  significance  are  the 
large  variations  in  apparent  signal  strength  which  can  be  caused  by  slight 
changes  in  the  aperture  to  object  separation.  Finally,  noise  can  be 
introduced  when  scanning  thick  OlOOOA),  translucent  samples,  because  the 
system  will  detect  light  scattered  from  the  diffraction  limited  far  field 
regions  as  well  as  from  the  collimated  near  field  regime.  In  short, 
considerable  effort  must  be  expended  to  increase  signal  strength  and  reduce 
noise  so  as  to  permit  rapid  high  resolution  microscopy. 

INITIAL  FEASIBILITY  STUDIES 

Our  initial  goal  wa>s  to  construct  a  simple  prototype  to  demonstrate  the 
feasibility  of  the  NSOM  concept.  The  essential  features  of  this  prototype 


are  discussed  below. 
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Subraicron  Aperture  Fabrication 


Aperture  Arrays 

Two  radically  different  techniques  have  been  employed  in  fabricating 
subraicron  apertures  for  near  field  microscopic  studies.  The  first  scheme 
consists  of  the  formation  of  aperture  arrays  in  a  thin,  planar  silicon  nitride 
membrane.  These  arrays  were  formed  using  electron  beam  lithography  at  the 
National  Research  and  Resource  Facility  for  Submicron  Structures.  This 
technique,  which  was  described  in  an  earlier  paper^7,  has  generated  uniform, 
reproducible  holes  as  small  as  80A  in  diameter.  After  aperture  information, 
the  opacity  of  the  surrounding  membrane  was  increased  by  evaporating  aluminum 
over  the  arrays.  Similar  apertures  have  already  been  applied  in  optical 
transmission  studies  in  which  we  have  demonstrated  that  visible  light  from  a 
tungsten  source  is  transmitted  through  holes  as  small  as  300A.3  The 
efficiency  of  the  transmission  process  is  demonstrated  in  Figure  3A,  where 
light  transmission  is  shown  through  2400A  diameter  apertures.  Figure  3B 
demonstrates  the  alteration  in  the  signal  to  noise  ratio  in  going  to  a  600A 
diameter  aperture  array.  The  uniformity  of  the  transmitted  intensity  attests 
to  reproducibility  of  this  method  of  aperture  fabrication.  In  addition  to  the 
utility  of  such  arrays  in  characterizing  the  optical  properties  of  submicron 
apertures,  they  can  be  applied  to  non-scanning  near  field  imaging  methods. 


Pipette  Apertures 

The  principal  problem  in  near  field  imaging  is  the  finite  depth  of  field 
» 

due  to  the  limited  extent  of  the  colliraation  of  radiation  in  the  near  field. 
Because  of  this  limitation,  apertures  in  a  planar  membrane  cannot  probe  reces¬ 
sed  regions  in  rough  surfaces.  To  address  this  problem,  we  have  developed  a 
simple,  inexpensive  method  of  aperture  fabrication  involving  the  formation  of 
subraicron  apertures  in  the  tips  of  highly  tapered,  metallized  glass  pipettes. 
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This  method  is  based  on  an  extension  of  the  technology  developed  for  making 
glass  pipettes  used  in  electrically  probing  cells  by  the  patch  clamp  method.18 
We  have  extended  this  technology  by:  a)  fabricating  apertures  with  diameters 
of  less  than  1000A,  b)  characterizing  the  apertures  at  the  tips  of  the 
pipettes  with  scanning  electron  microscopy,  c)  demonstrating  that  light  is 
transmitted  through  such  pipette  apertures. 

A  two  stage  method  was  used  for  constructing  the  apertures.  The  first 
stage  tapered  the  pipette  and  the  second  stage  pulled  until  the  breakage  point 
was  reached.  By  varying  the  wall  thickness  of  the  glass  pipettes,  aperture 
diameters  from  less  than  100QA  to  5000A  with  outer  diameters  of  5000A  to  7500A 
were  reproducibly  generated.  Thicker  walled  pipettes  were  found  to  produce 
smaller  apertures.  Aluminum  was  again  evaporated  to  increase  the  opacity  of 
the  glass  wall  surrounding  the  aperture.  This  evaporation  resulted  in  a 
uniformly  metallized  tapered  pipette  with  an  aperture  smaller  than  1000A  at 
the  tip  (see  Figure  4). 

The  taper  of  the  pipette  is  of  particular  importance.  When  light  is 
passed  down  the  pipette,  it  is  transmitted  through  the  outer  glass  walls  as 
well  as  through  the  central  region  (see  Figure  5).  The  inner  diameter  of  the 
pipette  rapidly  tapers  to  dimensions  of  less  than  an  optical  wavelength.  This 
model  of  light  transmission  suggests  that  the  radiation  exists  in  a  propagat¬ 
ing  mode  throughout  the  length  of  the  pipette,  since  the  pipette  outer  dia¬ 
meter  is  larger  than  the  cut-off  value  calculated  by  treating  it  as  a  classi¬ 
cal  wave  guide.  Only  at  the  thin  metallized  region  at  the  tip  of  the  pipette 

/ 

is  the  cut-off  threshold  reached.  Thus,  the  large  throughputs  observed  are  a 
consequence  of  the  fact  that  the  region  of  evanescence  is  short.  Upon  leaving 
the  aperture,  the  radiation  once  again  exists  m  a  propagating  mode. 
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An  NSOM  Prototype — Design 

The  pipette  type  apertures  were  an  important  advance  in  our  development 

of  a  simple  NSOM  prototype,  but  even  for  the  initial  feasibility  experiments, 

the  technical  challenges  discussed  above  had  to  be  addressed  in  at  least  a 
preliminary  fashion.  The  consideration  of  these  challenges  in  our  development 
of  a  workable  NSOM  prototype  will  now  be  discussed. 

Initial  System 

A  two-stage  system  was  employed  to  shield  the  instrument  from  floor 

vibrations.  The  first  consisted  of  a  large  table  top  mounted  on  tennis 

balls.  Measurements  in  our  laboratory  demonstrated  that  this  system  reduced 

high  frequency  vertical  disturbances  to  amplitudes  of  less  than  300A.  The 

second  consisted  of  a  heavy  steel  plate  resting  on  a  foam  pad,  which  provided 

isolation  in  the  intermediate  frequency  range  (10  to  30  Hz). 

Using  precision  translation  stages  and  conventional  micrometers,  the 

pipette  was  positioned  with  approximately  one  micron  accuracy  in  the  xy  plane 

over  the  section  of  the  sample  to  be  scanned.  Another  stage  was  used  to 

position  the  pipette  within  a  few  microns  of  the  sample  surface.  To  bring  the 

pipette  within  the  near  field  region,  a  stacked  piezoelectric  element  was  used 

(see  Figure  6).  A  low-noise  (15raV  peak  to  peak),  high-voltage  supply  was 

employed  to  insure  that  the  error  in  the  piezoelectric  expansion  due  to 

voltage  ripple  was  considerably  less  than  the  20A  limit  mentioned  previously. 

Another  piezoelectric  element  was  used  to  drive  the  pipette  in  the  scan 

direction.  The  step  size  could  be  varied,  but  was  typically  150A.  Because 

/ 

conventional  piezoelectric  elements  have  limited  expansion  ranges,  the  scans 
were  constrained  to  be  no  larger  than  several  microns.  Furthermore,  these 
conventional  elements  exhibit  hysteresis  problems  (=152  of  the  full 
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expansion),  so  that  accurate  positioning  could  be  achieved  only  along  a  single 
axis. 

To  deal  with  the  problems  posed  by  signal  to  noise  considerations,  the 
prototype  system  (see  Figure  6)  employed  an  optical  multichannel  analyzer  with 
an  SIT  vidicon  detection  system  in  conjunction  with  a  conventional  light 
microscope.  This  design  insured  that  a  sufficient  portion  of  the  available 
light  was  collected  and  detected.  Noise  in  each  channel  of  the  detector  was 
at  most  half  the  signal  due  to  the  aperture,  even  when  small  apertures  and 
weak  (tungsten)  illumination  were  used.  Much  of  the  noise  consisted  of  a 
constant  background  level,  which  was  easily  subtracted.  To  measure  the  noise 
from  other  sources,  the  final  detection  signal  was  plotted  as  a  function  of 
time  in  a  test  situation  where  the  aperture  position  was  fixed.  The  resulting 
variations  in  intensity  were  roughly  one  percent  of  the  apparent  signal. 

These  noise  fluctuations  were  at  least  an  order  of  magnitude  smaller  than  the 
intensity  changes  observed  in  the  course  of  a  super-resolution  scan. 

NSOM  Prototype — Results 

Microfabricated  Patterns  for  Resolution  Tests 

To  test  the  NSOM  method,  it  was  necessary  to  fabricate  objects  that  could 
be  used  to  accurately  define  the  resolution  limits  of  the  resulting  device. 
Toward  this  end,  electron  beam  lithography  was  used  to  create  a  mask 
consisting  of  several  sets  of  line  gratings  with  varying  periodicity.  The 
image  of  this  mask  was  reproduced  with  high  fidelity  as  a  series  of  aluminum 
lines  on  a  glass  coverslip  using  the  method  of  contact  printing.  These  lines 
were  2000A  wide  and  were  separated  by  spacings  ranging  from  2000A-8000A.  In 
Figure  7,  a  scanning  electron  micrograph  of  one  of  these  gratings  is  shown. 
Note  the  quality  with  which  the  pattern  is  reproduced.  As  noted  above  such 
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high-quality  test  objects  with  reproducible  features  are  required  before  any 
conclusions  can  be  drawn  about  NSOM  resolution. 

A  One  Dimensional  Scan 

Using  the  NSOM  prototype  and  the  test  object  discussed  above,  a  one 
dimensional  scan  was  performed  which  demonstrated  super-resolution 
capabilities.  The  light  from  a  100W  tungsten  lamp  was  passed  through  a 
pipette  to  illuminate  the  aperture.  The  tip  was  brought  into  contact  with  the 
sample  surface  and  then  retracted  several  nanometers  for  scanning. 

The  initial  scan  was  performed  using  a  5000A  diameter  pipette  type 
aperture  and  the  results  are  shown  in  Figure  8.  Note  that  the  periodicity  of 
the  gracing  pictured  in  Figure  7  is  reproduced  throughout  the  entire  length  of 
the  NSOM  scan.  A  resolution  of  better  than  1500A  is  inferred  from  averaging 
the  sharpness  of  the  steps  in  going  from  transmission  minima  to  adjacent 
transmission  maxima  and  vice  versa.  Note  further  that  the  technique  permits 
detection  of  sharp  edges  with  a  resolution  considerably  better  than  the  limit 
expected  based  on  a  naive  model  of  the  aperture  size.  This  is  understandable, 
because  with  a  sufficient  signal  to  noise  ratio  it  should  be  possible  to 
detect  changes  in  intensity  due  to  partial  occlusion  of  the  aperture  by  the 
sharp  edge. 

Future  Improvements 

Although  the  prototype  system  was  sufficient  to  demonstrate  the 
super-resolution  concept,  it  suffered  from  a  number  of  limitations.  These 
included  the  inability  to  3can  in  two  dimensions,  the  inability  to  scan  over 
rough  topography,  the  inability  to  produce  highly  reproducible  scans,  and  the 
inability  to  scan  at  high  speed.  Methods  exist  for  attacking  all  these 
problems.  We  are  currently  incorporating  these  methods  in  the  construction 
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of  a  second  generation  NSOM  instrument.  Many  of  the  ideas  employed  in  this 
construction  effort  were  first  developed  for  scanning  tunneling  microscopy, 
where  positioning  on  a  sub-Angstrora  scale  is  required.  At  the  completion  of 
this  effort,  we  hope  to  have  a  versatile  system  capable  of  investigating 
problems  of  biophysical  interest. 

BIOPHYSICAL  APPLICATION 
Fluorescence  Detection 

Because  NSOM  employs  visible  radiation,  many  of  the  techniques  which  have 
already  been  developed  for  conventional  microscopy  (e.g.,  polarized 
illumination,  phase  contrast  methods,  etc.)  can  be  used  to  advance  the 
biophysical  applicability  of  super-resolution  microscopy.  In  particular, 
contrast  enhancement  by  fluorescence  techniques  will  be  useful,  because  many 
biological  systems  can  be  selectively  labeled  with  fluorophores  in  a 
non-destructive  fashion  or  else  contain  fluorescent  domains  in  their  natural 
state.  The  only  potential  problem  associated  with  fluorescent  labeling  lies 
in  the  reduction  of  the  signal  strength  due  to  the  limited  cross  sections  and 
quantum  yields  of  these  fluorophores.  However,  in  this  paper,  we  have  taken  a 
first  step  in  demonstrating  the  feasibility  of  fluorescence  imaging  by 
observing  fluorescence  emanating  from  sub-wavelength  apertures  (see  figure 
9).  This  result  was  obtained  by  illuminating  a  thin  film  of  perylene  with 
5145A  radiation  from  a  laser,  and  viewing  the  fluorescence  transmitted  through 
an  aperture  array  with  the  appropriate  filter.  Note  that  relatively  strong 
signals  are  obtained,  ev,en  though  perylene  has  only  a  moderate  quantum 
efficiency. 

Fluorescence  From  a  Functioning  Biological  System 


The  next  step  in  approaching  a  significant  biological  problem  involved 
observing  the  transmission  of  fluorescent  light  through  subraicron  apertures 
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eraenating  from  a  viable  biological  specimen.  To  acheive  this,  spinach 
chloroplaats  were  used.  A  suspension  of  such  chloroplasts  was  placed  directly 
on  an  array  of  2400A  holes,  and  the  suspension  was  made  sufficiently  dense  so 
that  the  chloroplasts  adhered  to  the  array  and  did  not  migrate.  The 
chloroplasts  were  illuminated  through  the  apertures  using  4579A  laser 
radiation  incident  from  the  opposite  side  of  the  array,  and  the  fluorescence 
was  detected  back  through  these  same  apertures  in  reflection.  The  results, 
shown  in  Figure  10,  demonstrate  the  feasibility  of  detecting  such  fluorescence 
on  a  subraicron  scale  in  a  non-destructive  manner. 

Chloroplasts  were  chosen  for  out  initial  experiments  because  they  are 
involved  in  the  current  debate  concerning  the  distribution  of  the  two 
photosysteras  involved  in  photosynthesis.  Photosystem  II  is  involved  with  the 
splitting  of  water  while  photosystem  I  is  involved  with  the  alteration  in  the 
oxidation-reduction  potential  of  electrons  ejected  from  the  split  water.  Both 
photosystems  are  found  in  the  thylakoid  membranes  of  chloroplasts,  but  only 
photosystera  II  fluoresces  at  room  temperature. 

The  membranes  themselves  contain  two  structurally  distinct  regions:  the 

grana,  in  which  the  membranes  are  stacked;  and  the  stroma,  which  are  regions 

of  unstacked  membranes  interconnecting  the  grana  regimes.  According  to  one 

model,  photosystera  I  is  excluded  from  the  stacked  regions  of  the  grana. 

However,  the  evidence  on  which  this  model  is  based  comes  from  freeze-fracture 

electron  microscopy,  detergent  methods  of  fractionating  the  grana  and  stroma 
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membranes,  and  mechanical  fractionating  methods  1  ,  all  of  which  have  a 

/ 

disruptive  effect  on  the  membranes.  Conventional  microscopy  cannot  be  used, 
since  the  grana  regions  cannot  be  resolved  with  sufficient  clarity  to 
determine  the  distribution  of  the  two  photosysteras.  In  contrast,  near  field 
microscopy  should  provide  the  resolution  necessary  to  determine  the 
photosystera  distribution  within  viable  chloroplasts  maintained  in  their  native 


environment.  We  are  currently  working  on  this  problem,  the  results  of  which 
should  help  elucidate  energy  transfer  mechanisms  involved  in  photosynthesis. 
Signal  to  Noise  Ratio  Calculations  for  Biological  Applications 

The  signals  detected  by  fluorescence  NSOM  imaging  are  weak  and  many 
possible  sources  of  noise  can  be  expected.  In  this  section  we  demonstrate  the 
feasibility  of  fluorescence  imaging  methods  by  calculating  the  signal  to  noise 
ratio  in  a  specific  biological  application.  The  particular  problem  we  focus 
on  is  the  detection  of  a  receptor  in  Che  plasma  membrane  of  human  fibroblast 
cells  that  binds  low  density  lipoprotein  (LDL).  LDL  is  the  predominant 

transporter  of  plasma  cholesterol,  and  at  least  30  fluorescent  lipid  molecules 
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dan  be  incorporated  per  LDL  molecule.  It  has  been  shown  that  with  the 
fluorescent  lipid  3,  3 '-dioctadecylindocarbocyanine  iodide  (dil[3])  single  LDL 
molecules  on  single  LDL  receptors  can  be  detected.  We  can  use  this  data  to 
obtain  signal  to  noise  estimates  which  will  indicate  the  feasibility  of 
imaging  such  receptors.  However,  these  calculations  are  not  limited  to  this 

example,  because  phycobiliproteins  have  been  shown  to  be  as  highly  fluorescent 

2  3 

as  LDL  and  can  be  used  to  selectively  label  a  variety  of  cellular  receptors. 

We  have  been  able  to  irradiate  LDL  molecules  with  green  light  for  about  one 

minute  at  an  intensity  of  lpW/p  before  the  resulting  fluorescent  image  begins 

22  2 

to  fade.  This  corresponds  to  a  total  dose  of  1.9x10  photons/cra  . 

To  determine  how  this  limit  affects  the  incident  intensity  and  hence  the 
total  signal  in  an  NSOM  application,  we  consider  here  the  specific  case  of 
near  field  microscopy  in  reflection.  The  use  of  such  a  geometry  is 
advantageous,  because  e&ch  molecule  is  affected  only  when  the  illuminated 
aperture  is  overhead.  In  practice,  this  is  not  strictly  true,  because  the 
finite  opacity  of  the  screen  must  be  taken  into  account.  For  example,  a  600A 
thick  layer  of  aluminum  will  transmit  roughly  0.01Z  of  the  incident  light.  If 
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the  attenuation  of  the  radiation  in  the  aperture  is  assumed  to  be  90%,  then 
the  differential  transraittability  of  the  aperture  relative  to  the  screen  is 
103  . 

In  a  typical  scan,  we  may  wish  to  record  a  100x100  element  image  within 

100  seconds.  In  this  case,  the  differential  transraittability  indicates  that 

the  LDL  molecules  will  receive  dosages  roughly  ten  times  greater  from 

transmission  through  the  screen  as  they  will  from  transmission  through  the 

aperture.  Hence,  to  remain  below  the  photobleaching  threshold  within  the  100 

second  scan  time,  a  maximum  photon  flux  past  the  screen  of 
20  2 

1.7x10  photon/sec/cm  is  permissible,  corresponding  to  an  incident  intensity 
of  5.4raW/ii2. 

Although  at  first  glance  it  may  appear  troubling  that  the  dosage  from 

screen  transmission  is  an  order  of  magnitude  greater  than  that  from  aperture 

transmission,  when  the  situation  is  considered  in  detail  it  is  not 

problematical  at  all.  First,  it  must  be  remembered  that  these  figures  refer 

to  the  total  dosages  in  the  course  of  a  scan,  and  at  any  given  time  Che 

intensity  directly  below  the  aperture  is  roughly  1000  times  greater  than  the 

intensity  below  the  screen.  Furthermore,  from  the  standpoint  of  the  signal  to 

noise  ratio,  we  are  only  interested  in  the  fluorescent  light  which  is 

transmitted  back  through  the  aperture  and  screen  for  eventual  detection. 

In  as  much  as  this  is  the  case,  the  transmission  of  light  through  the  screen 

is  unimportant,  because  any  fluorescence  which  might  otherwise  be  generated  by 

LDL  raolecul  s  under  the  screen  is  highly  quenched  by  the  presence  of  the  metal 
24  ' 

within  the  screen.  Hence,  only  the  fluorescence  from  Che  near  field 
directly  below  the  aperture  contributes  to  the  signal,  and  only  far  field 
fluorescence  contributes  to  the  noise. 

The  maximum  permissible  intensity  calculated  above  can  now  be  used  to 
estimate  the  total  signs,  which  can  be  expected  in  our  Gedanken  experiment. 


In  particular,  by  again  assuming  a  10Z  transraittabi lity  for  the  aperture  and 

noting  that  the  above  scan  parameters  lead  to  a  10ms  exposure  time  per  pixel 
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location,  we  find  that  1.7x10  photons/cm  for  each  pixel  are  transmitted 
through  the  aperture.  We  will  further  assume  that  these  photons  are  all 
incident  on  the  fluorescent  site  of  a  single  receptor. 

To  calculate  the  resulting  signal,  we  note  that  each  of  the  30  dil 
chromophores  per  LDL  molecule  has  an  effective  photoabsorption  cross  section 
of  3.8x10"  cm  with  a  fluorescent  yield  of  0.1.  Using  these  numbers,  we 
find  an  initial  signal  of  1.9x10^  fluorescent  photons/pixel. 

A  number  of  factors  cause  the  detected  signal  to  be  lower  than  this 
initial  value.  First,  in  the  reflection  mode  described  here,  the  assumed 
aperture  transmission  reduces  the  fluorescent  light  intensity  by  90Z.  Second, 
inefficiencies  in  the  collection  optics  can  be  assumed  to  further  reduce  the 
intensity  by  50Z.  Finally,  the  photomultiplier  that  will  form  the  heart  of 
a  sensitive  detection  system  has  a  quantum  efficiency  of  only  20Z.  When  all 
these  factors  are  taken  into  account,  the  detected  signal  is  found  to  be 
i.9xl03  counts/pixel. 

To  determine  the  noise  associated  with  a  signal  of  this  strength,  we 
first  consider  noise  due  to  photon  counting.  Statistical  noise  for  the  above 
counting  level  is  roughly  44  counts/pixel.  In  contrast,  dark  noise  from  the 
photomultiplier  should  be  less  than  10  counts/sec.  For  the  above  scan  rates, 
this  corresponds  to  0.1  counts/pixel,  and  hence  is  quite  negligible. 

Other  fundamental  sources  of  noise  arise  from  the  limited  extent  of  the 
collimation  of  radiation  in  the  near  field.  For  example,  although  the  LDL 
molecules  in  close  proximity  to  the  metal  screen  will  not  exhibit 
fluorescence,  other  such  molecules  bound  to  the  far  side  of  the  cell  membrane 
will  be  weakly  illuminated  by  residual  far  field  radiation  from  the  aperture, 
and  hence  may  yield  dim  fluorescence  noise.  This  noise  is  considerably  less 
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than  the  signal  due  to  fluorescence  from  receptors  within  the  intense  near 
field  radiation  pattern  emanating  from  the  aperture,  as  can  be  understood  by 
the  following  simple  argument. 

As  was  mentioned  previously,  the  intensity  of  radiation  emanating  from  a 
slit  falls  in  an  almost  exponential  manner  within  the  near  field.  It  is 
reasonable  to  assume  that  a  similar  behavior  exists  in  the  case  of  a  circular 
aperture.  In  fact,  by  treating  an  aperture  of  diameter  d=500A  as  a  waveguide, 
a  specific  form  of  exp[-2nr/d]  can  be  assumed  to  hold  out  to  a  colliraation 
distance  of  roughly  250A  .  Hence,  the  radiation  should  remain  collimated 
through  the  entire  thickness  of  that  portion  of  the  cell  membrane  which  is 
directly  under  the  aperture,  and  any  fluorescence  that  occurs  in  this  region 
can  be  considered  to  contribute  to  the  signal.  In  contrast,  any  fluorescence 
which  occurs  in  the  far  field  region  contributes  to  the  noise.  However, 
within  the  far  field,  the  aperture  can  be  modeled  as  a  dipole  radiator,  so 
that  the  outward  flux  falls  as  1/r6  .  If  we  assume  that  this  behavior  is  valid 
up  to  r*250A,  then  the  continuity  of  intensity  at  this  distance  can  be  used  to 
provide  a  complete,  qualitiative  model  of  the  energy  flow  past  the  aperture. 

We  can  now  use  such  a  model  to  estimate  the  signal-to-noise  ratio  due  to 

the  limited  extent  of  the  collimated  near  field.  The  signal  is  determined  by 

integrating  the  exponential  intensity  factor  over  a  hemisphere  of  250A 

radius.  Similarly,  by  integrating  the  factor  of  1/r6  over  a  250A  thick  shell 

one  micron  away  (which  represents  the  membrane  furthest  from  the  aperture), 

the  noise  is  found.  In  this  example,  which  should  fairly  well  model  the 

» 

binding  of  LDL  receptors  to  fibroblast  cells,  the  noise  is  found  to  be  several 
orders  of  magnitude  less  than  the  signal.  Even  if  we  were  to  assume  that 
chroraophores  are  evenly  distributed  throughout  the  bulk  of  the  cells  (which  is 
certainly  not  the  case  in  our  example),  a  signal  to  noise  ratio  of  7.7  can  be 
obtained  using  this  model.  When  it  is  also  considered  that  the  fluorescence 


from  the  chroraophores  decreases  in  intensity  by  1/r  as  well,  the  noise  is 
dropped  even  further. 

The  implications  of  this  last  point  cannot  be  overemphasized.  The  finite 
extent  of  the  near  field  is  used  to  advantage  since  the  far  field  noise  is 
small  compared  to  the  near  field  signal.  Thus,  NSOM  can  be  used  to  image 
surfaces  on  thick  samples  without  sectioning. 

Returning  to  our  model  for  imaging  LDL  receptors,  it  can  now  be  seen  that 
the  most  significant  source  of  noise  is  statistical  in  nature.  Thus,  by  using 
the  above  results,  we  find  that  a  signal-to-noise  ratio  of  about  44  should  be 
obtainable,  and  hence  the  experiment  should  be  entirely  feasible.  Indeed,  it 
would  not  be  necessary  to  run  the  experiment  at  the  threshold  for 
photobleaching.  If,  for  example,  the  incident  intensity  were  to  be  decreased' 
by  a  factor  of  ten,  a  signal  to  noise  ratio  of  almost  14  would  still  be 
achieved.  Alternatively,  the  size  of  the  scan  speed  could  be  increased  by  a 
factor  of  ten  with  an  identical  change  in  the  signal-to-noise  ratio. 

CONCLUSIONS 

Our  initial  prototype  system  has  demonstrated  the  feasibility  of  NSOM  as 

a  high  resolution  imaging  technique.  Furthermore,  the  results  of  the 

fluorescence  experiments  and  the  signal-to-noise  ratio  calculations  indicate 

that  NSOM  may  be  applied  to  a  variety  of  biological  problems.  For  example, 

with  the  development  of  a  flexible  second  generation  system,  NSOM  could  be 

25 

used  to:  image  lipid  microdomains  in  model  membrane  assemblies  or  living 
cells;  investigate  banding  patterns  in  chromosomes  with  fluorescent  molecules 

S 
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that  bind  to  particular  DNA  sequences  and  structures  ;  delineate 
cystoskeletal  interactions  with  membrane  bound  proteins;  characterize 

27 

alterations  in  cell  surface  topography  as  a  result  of  cellular  stimulation  ; 

28 

visualize  organelle  movement  along  dissociated  actin  filaments  ;  and  to  study 
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the  kinetics  of  specific  proteins  during  muscle  contraction.  In  addition, 
this  technique  could  be  used  to  investigate  problems  in  many  other  areas  from 
microelectronics  to  chemistry  where  non-destructive,  high  resolution  imaging 
is  required. 
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FIGURE  CAPTIONS 


Figure  1  .Schematic  representation  showing  the  colliraation  of  radiation 
emanating  from  a  sub-wavelength  aperture. 

Figure  2  The  energy  flux  emanating  from  a  500A  wide  slit,  normalized  to  the 
energy  of  the  incident  plane  wave,  is  shown  plotted  as  a  function 
of  both  the  perpendicular  distance  (z)  from  an  1800A  thick 
conducting  screen  and  the  ratio  of  the  lateral  distance  from  the 
center  of  the  slit  to  the  half  width  of  the  slit  (x/a).  Notice 
that  the  radiation  is  collimated  to  roughly  the  edge  of  the  slit 
(x/a*l)  out  to  approximately  z=,250A. 

Figure  3  (A)  5145A  light  transmission  through  2400A  diameter  apertures  is 

shown  detected  in  the  far  field  with  an  optical  multichannel 
analyzer.  Note  the  uniformity  of  the  light  transmission  through 
the  nine  apertures  shown.  Incident  beam  power  was  0.62ra  W/cra  . 

(B)  5145A  light  transmission  through  600A  diameter  apertures  is 

shown  as  detected  in  the  far  field  with  an  optical  multichannel 
analyzer.  Incident  laser  power  was  as  above.  Note  the  even  for 
these  smaller  apertures  a  good  signal-to-noise  ratio  is  achieved. 

Figure  4  A  scanning  electron  micrograph  of  the  tip  of  a  metallized  pipette 
is  shown  with  a  500QA  outer  diameter  and  with  a  <1000A  central 
aperture.  Although  even  smaller  apertures  could  be  prepared  with 
our  technique,  they  would  be  difficult  to  characterize  with  the 
limited  resolution  of  the  SEM  available. 

I 

Figure  5  A  cross-sectional  representation  of  a  metallized  pipette  is  shown 
together  with  the  propagation  of  light  toward  the  tip. 


23 


Figure  6 


Figure  7 


Figure  8 


Figure  9 


Figure  10 


A  schematic  diagram  of  the  experimental  apparatus  used  for 
our  initial  feasibility  experiments  is  shown.  Both  the  optical 
multichannel  analyzer  and  the  piezoelectric  translators  are 
controlled  with  a  microcomputer. 

A  scanning  electron  micrograph  of  a  grating  of  aluminum  lines  on 

glass  used  for  our  resolution  tests  is  shown.  The  lines  were  2000A 

separated  by  400 0A .  To  obtain  this  micrograph  the  sample  was 

coated  with  a  100A  thick  film  of  gold  palladium. 

A  scan  of  the  test  grating  in  Figure  7  is  shown.  From  the 

sharpness  of  the  intensity  modulation  between  adjacent  maxima  and 

minima,  a  resolution  of  <1500A  is  inferred.  Furthermore,  note  that 

the  periodicity  of  the  grating  is  reproduced  over  the  entire  scan. 

Perylene  fluorescence  through  (a)  1200A  diameter  and  (b)  600A 

diameter  apertues  are  shown.  The  exciting  laser  light  at  5145A 

illuminated  the  sample  from  below, and  the  fluorescence  was  detected 

through  the  apertures  from  above.  The  laser  power  that  was  used 

caused  no  detectable  bleaching  of  the  sample.  The  fluorescence  was 

recorded  with  ASA  400  Kodak  Ektachrome  film  using  a  30  sec 

exposure.  These  arrays  contain  graduated  aperture  sizes  with  the 

maximum  corresponding  to  the  diameters  noted  above. 

Photosystera  H  fluorescence  is  shown  detected  through  2400A 

diameter  apertures  using  epi-iliumination.  The  incident  intensity 

from  an  argon  ion  laser  at  4579A  was  0.067  dW/p  at  the  surface  of 
/ 

the  array.  Kodacolor  ASA  400  film  and  an  exposure  time  of  40  sec 
was  used  to  detect  the  fluorescence  from  the  chloroplasts . 
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